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ABSTRACT 

It is now well established that many galaxies have nuclear star clusters (NCs) 
whose total masses correlate with the velocity dispersion a of the galaxy spheroid in a 
very similar way to the well-known supermassive black hole (SMBH) M — a relation. 
Previous theoretical work suggested that both correlations can be explained by a 
momentum feedback argument. Observations further show that most known NCs have 
masses < 10 8 M Q , while SMBHs frequently have measured masses > 10 8 M Q , which 
remained unexplained in earlier treatments. We suggest here that this changeover 
reflects a competition between the SMBH and nuclear clusters in the feedback they 
produce. When one of the massive objects reaches its limiting M — a value, it drives 
the gas away and hence cuts off its own mass and also the mass of the "competitor" . 
The latter is then underweight with respect to the expected M — a mass. 

More specifically, we find that the bulge dynamical timescale is a steeply rising 
function of velocity dispersion, and that the NC-SMBH changeover occurs where the 
dynamical time is about equal to the Salpeter time. We propose that SMBHs, growing 
on the Salpeter time scale, arc unable to reach their M — a mass quickly enough in 
small bulges. The central regions of these bulges are swamped with gas which fragments 
into stars, creating the nuclear clusters. The latter then limit their own growth by the 
feedback they produce, settling on their (offset) M — a relation. The SMBH in such 
bulges should be underweight as their growth is curtailed before they reach the M — a 
mass. In large bulges, on the other hand, the SMBH catches up quickly enough to 
settle on its M — a relation. Nuclear star clusters may also exist in such bulges but 
they should be underweight with respect to their M — a sequence. 
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INTRODUCTION ter. For an isothermal potential, M{R) and M tot ai(-R) 

proportional to R, so the result is 



It is well known that the masses of the supermassive black 

holes (SMBHs) in the nuclei of early-type galaxies and yy _ ^fg a ^\ 
bulges correlate with the velocity dispersions of the stellar G 

spheroids (e.g.. iGebhardt et all 12000: |Ferrarese fc MerrittJ Here f g is the baryonic fraction and a 2 = GM/2R is the 

l200d: iTremaine et all 120021') A s i mple explanation invokes velocity dispersion in the bulge. To order of magnitude, the 

momentum feedback (|Kindl2003l. 120051). In this picture the relation [2] holds for any potential if estimated at the virial 

SMBH luminosity is limited by the Eddington value, and radius. 

the momentum outflow rate produced by radiation pressure Requiring that momentum output produced by the 

is of the order of black hole should jus t balance the weig ht of the gas leads to 



^Edd 4tyGMbh 



black nole should just balance the weigh 
the AfsH-o- relation (|Kindl2003l . l2005h : 



IISMBHW— — = , (1) f K 4 

C K Mbh = — T^CT . (6) 

7TG 

where k is the opacity, assumed to be dominated by the The model is attractive in its physical simplicity. Further, 

electron scattering, and Mbh is the SMBH mass. This mo- the result contains no free parameters, but is very close to 

mentum flux produces an outward force on the gas in the the observed Mbh-o" relation. 

bulge, whose weight is W(R) = GM(R)[M totill (R)]/ R 2 , Another feature commonly found in the centres of 

where M(R) is the enclosed gas mass at radius R, and galaxies are nuclear s tar clusters. They are found in 

Mtotal 



lU(li) is tne enclosed gas mass at radius li, and galaxies are nuclear s tar clusters, iney are tound m 
(R) is the total enclosed mass including dark mat- late type spirals (e.g., iBoker et al.l l2002l h bulgeless spi- 
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rals (|Walcher et al.ll2005l ). ed ge-on spirals dSeth et al.ll2006h , 
and dwarf elliptical galaxies ( Cote et al.ll2006l ). The cluster 
masses range from 10 6 to 10 s Mq (althoug h even more mas- 
sive clusters have been found recently by iKormendv et all 
l2009h . whereas their sizes are only a few to a few tens parsec. 
Intriguingly, for dwarf ellipticals the masses of NCs corre- 
late with the host properties. Namely, iFerrarese et al.l <|2006h 
found that the NC mass is related to the bulge velocity dis- 
persion in exactly same way as the SMBH - velocity dis- 
persion relation, but with normalisation offset by about an 
order of magnitude. The cluster masses are also almost lin- 
early proportional to the tot al bulge mass (|Ferrarese et al.l 
120061; IWehner fc Harrij|2006h . 

iMcLaughlin et aL ( 20061 ) proposed that the observed 
Mnc — o relation for dwarf elliptical galaxies f ollows natu- 
rally from an extension of the above arg ument |Kind |2003l . 
120051 )) to the outflows from young star clusters containing 
massive stars. These individual stars are also Eddington- 
limited, and produce outflows with momentum outflow rate 
~ -^Edd/c where LEdd is calculated from the star's mass. 
Young star clusters with normal IMFs produce momentum 
outflow rate 

n NC « (4) 

c 

where A w 0.05 and LEdd is now formally the Eddington 
value corresponding to the total cluster mass. To produce 
the same amount of momentum feedback, a young star clus- 
ter must therefore be 1/A times more massive than a SMBH 
radiating at the Eddington limit, and hence: 

/fl« 4 



AL 



XttG 2 



(5) 



Strikingly, 1/ A is quite close to the offset in mass between the 
A^bh-o" and Mnc -0 " rel ations. Furthermore, w hile equation 
[T]is very plausible fsee lKing fc Pounds! 120031 ). the momen- 
tum outflow rates from stars are known in detail observa- 
tionally. The largest uncertainty in the equation [3] is there- 
fore t he stellar IMF , whi ch is observationally fairly constant 
fe.g.. lKroupall2002l ). The IMcLaughlin et alj (|200ei ) explana- 
tion of Mnc-o - re lation thus appears similarly robust to the 
iKind |2003l . 120051) model for SMBH fee dback. 



However, 



McLaughlin et al.l (|2006l ) did not offer an ex- 



planation of why bulges with smaller a contain nuclear clus- 
ters, while more massive galaxies contain SMBHs and not 
NCs. Here we propose an explanation, noting that timescales 
are important in this problem as well as energetics. 

Our simple theory for the observed bimodality of NC 
and SMBHs is based on the premise that the dominant ob- 
ject must be able to grow quickly and yet stay active for 
long enough to provide the needed feedback. As we show 
below, in small bulges this argument favours nuclear star 
clusters whereas in larger ones the situation is reversed. Be- 
low we explain our idea, address observational constraints 
and conditions needed for it to work, and suggest possible 
astrophysical implications. 

We note that nuclear star clusters could in principle 
form elsewhere in the galaxies and then migrate inwards due 
to dynamical frictio n with the background stars. However, 
iMilosavlievidl |2004 ) argues against this possibility due to 
the short time scales available for this process, and argues in- 
stead that these clusters may form in situ. We agree with this 
point, and further note that observations of young massive 



stars in the central parsec of the Milk y Way offer direct sup- 
port t o the in-situ formatio n model jNavakshin fc Cuadral 
120051 ; iPaumard et al. I 120061 : iNavakshin fc Sunvaevl 120051 ). 
The exact geometrical arrangement of the forming stars (a 
thick disc or a quasi-spherical cluster) is irrelevant on the 
scales of the parent galaxy. 



2 TIMESCALES 

Black holes and nuclear clusters each evolve on character- 
istic timescales. SMBH growth is limited by the Eddington 
accretion rate, MEdd = LEdd/(ec 2 ), where e ~ 0.1 is the 
radiative efficiency of accretion. SMBH masses can grow no 
faster than exp(f/isai P ), where 

AfBH Kec 

= 



4.5 x 10 £0.1 yr 



(6) 



M Ed d 

is the Salpeter time, with eo.i = e/0.1. Star formation can 
occur on the free-fall or dynamical timescale tdyn of the 
system, which is less tha n a million years fo r many observed 
young star clusters (e.g.. lHillenbrandlll997f l. 

Once a SMBH is created, its feedback can be activated 
at any time, provided that the accretion rate is high enough. 
By contrast, star cluster feedback has a "half life" of around 
tMS ~2 x 10 7 yr, since this is the main-sequence lifetime 
of the massive stars contributing most to the feedback (e.g. 
iLeitherer et al.lll992f l. This timescale is only a factor of two 
shorter than the Salpeter time. After this time (t = i m a), the 
ability of the nuclear clusters to expel gas from the galaxy 
is severely reduced. They would have to be rebuild their 
population of massive young stars to restart. It is not obvi- 
ous that this is physically possible inside an existing dense 
stellar cluster. 

Consider a bulge where the dynamical time (equation [7] 
below) is much shorter than the Salpeter time. In a gas feed- 
ing event (e.g. a merger) , the bulge regains dynamical equi- 
librium before any significant SMBH growth and feedback 
sets in, and there is nothing to prevent gas from collecting in 
the bulge centre. The accumulated gas is then consumed by 
star formation in nuclear regions, forming nuclear clusters 
which quickly reach their Mnc _ o" limiting mass. This cuts 
off growth of everything - the bulge, the NC and the SMBH 
as well. The SMBH in these bulges are thus bound to be 
underweight compared to the Mbh-c" relation. 

In the opposite extreme, when the bulge dynamical time 
is longer than the Salpeter time, the SMBH can grow quickly 
enough to reach its limiting Mbh-ct mass. While nuclear 
star clusters might be created there as well, their feedback 
quickly (i.e. in about 20 million years) becomes negligible. 
The situation is thus the reverse of the last paragraph, and 
it is the nuclear star clusters that are underweight in these 
bulges. 

Below we estimate td yn as a function of bulge mass or ve- 
locity dispersion. We find that ts a i P > £dyn in smaller bulges 
(a < 150 kms -1 ) and tsaip ;$ tdyn in larger ones. 



DYNAMICAL TIME AND VELOCITY 
DISPERSION 



A tac i t but obv ious assumption in the a rguments of IKind 
l|2003l . l2005l ) and lMcLaughlin et al.l (|2006l ) is that the source 
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of feedback can respond quickly to bulge growth and thus 
influence it. Depletion of gas in the bulge by star formation 
and the onset of stellar feedback probably occur within a 
few bulge dynamical time scales 



t - * 



(7) 



where R and a are the scale length and velocity dispersion 
of the bulge. If the feedback source fails to reach its limiting 
mass its feedback remains unimportant. 

We now consider how the dynamical time td yn scales 
with the a of the stellar component of a galaxy. It 
is well established that stellar spheroids occupy a two- 
dimensional 'fundamental plane' in the space defined by the 
total luminosity, the scale length and the velocity d isper- 
sion (|Diorgovski fc Davisl Il9871 ; iBernardi et alj l2003h . The 
plane is tilted relative to its position expected from the 
virial theorem for a homologous population of galaxies in 
a dynamical equilibrium. Several explanations for the ori- 
gin of the tilt have been discussed in the literature, in- 
cluding systematic variations in the mass to light ratio of 
the stellar populations, or ch anges in the dark matter frac- 
tion (e.g. iDekel fc Coxl 120061 ). Recent studies have shown 
that the tilt is essentially independent of the wavelength of 
the observations, suggesting that stellar population varia- 
tions are not the dominan t contribution (|La Barbera et all 
120081 ; IBernardi et al.ll2003h . 

Projections of the fundamental plane lead to a num- 
ber of simple scaling relations. Using Sloan Digital Sky Sur- 
vey photometric and spectroscopic data for ~ 9000 galax- 
ies with measu r ed ve locity dispersions of 100-400 kms - , 
IBernardi et al.l (|2003l ) derived the following relations be- 
tween R, a and the total luminosity L: 



R = 2.6 



150 



1.6 x 10 10 L 



0.7()4±(].()2.-, 



kpc 



L 



1.6 x 10 10 L 



0.23±0.012 



kms 1 



(8) 
(9) 



(The power- law indices are taken from IDekel fc Coxl (|2006l ) 
who report updated values obtained by Bernardi (priv. 
comm.) using revised SDSS photometry.) These relations 
lead to another one between the dynamical time and the 
total luminosity, given by 



tdyn = — = 17 



L 



1.6 x 10 10 L 



Myr 



(10) 



Combining relations © and (|10p . we obtain the dynamical 
time as a function of velocity dispersion, namely 



idyn — 17 



150 kms 



Myr 



(11) 



We have taken a ~ 150 km s -1 as our fiducial value as 
the observations show that no nuclear clusters have been 
observed in systems with a > 150 kms -1 . The above re- 
lations show that this roughly coincides with the transi- 
tion between systems with dynamical times longer than the 
Salpeter time ((©J . The details of the transition may depend 
on the merger history of the galaxy. 



4 DISCUSSION 

We have seen that momentum feedback gives a simple phys- 
ical explanation of why galaxy bulges are dominated by nu- 
clear clusters for low velocity dispersions and by supermas- 
sive black holes for high dispersions. We have emphasised 
that given an injection of gas, e.g. from a merger, galaxies 
with dynamical times shorter than the Salpeter time cannot 
grow their central black holes sufficiently quickly to affect 
the gas infall. Gas accumulating in the central regions can- 
not cool and condense indefinitely, so nuclear star clusters 
form and produce feedback. The masses of these clusters 
saturate at the mass ((5]) when they expel the remaining gas. 
The hole thus remains close to its 'seed' mass, which is pre- 
sumably less than the value ©. 

Note that this line of argument does not imply that 
dwarf elliptical galaxies with low velocity dispersion do not 
contain massive black holes. We only suggest here that their 
growth is slow; it is not entirely forbidden. Therefore, these 
galaxies may still contain underweight SMBH, i.e., black 
holes with mass significantly less than the corresponding 
Mbh-ct value. It seems difficult to avoid building up a mas- 
sive black hole in the very centre of the cluster and galaxy 
potential well. Mergers of low mass holes may provide an 
inter esting additional window for gravitational wave astron - 
omy (|Matsubavashi et ai]|2004l ; [Amaro-Seoane et al"1l2009l '). 

There is an important constraint for our model to be 
applicable. Star formation in the inner parts of galaxies 
with higher velocity dispersions is not ruled out by the 
considerations of this paper. These galaxies could thus po- 
tentially build up nuclear clusters in their centres. If this 
happens faster than the Salpeter time, then the masses of 
these clusters should saturate at the value (|3]). Presumably, 
if the SMBH continues to grow at the Eddington rate, it 
could then reach its limiting Mbh~& mass. However this 
picture would predict very massive nuclear star clusters (up 
to Mnc ~ 10 10 Mq for a ~ 300 km/sec) which are not ob- 
served. 

This suggests that SMBH growth by accretion should 
be the dominant process in the central parsecs of galax- 
ies, and that star formation occurs only as an alternative 
when gas cannot be consumed by the hole quickly enough. 
The latter naturally occurs if (a) the material is first de- 
posited into the disc on small scales where star formation 
does not occur due to the strong SMBH tidal eff ect (e.g., 
iKolvkhalov fc Sunvaev|[l98b1 ; lKing fc Pringlell2007ft . and (b) 
the hole is fed at a super-Eddington rate. Then the hole ac- 
cretes the gas at the Eddington rate, expelling the rest. A 
good fraction of the expelled gas would probably not travel 
very far from the centre of the galaxy. Gas is likely to be 
expelled with a range of velocities, some too low to escape 
to infinity. As its angular momentum is very low, the gas 
can fall back into the inner parsec(s). Such effects are ac- 
tu ally observed in the simulations of accretion discs winds 
bv lProgal (|2003l ). Deposited back into the accretion disc on 
parse c scales, the gas would then fuel star formation there 
(e.g., lGoodmarJl2003l ). 

Summarising this, applicability of our model demands 
that SMBH feeding be primary and star formation sec- 
ondary in the inner few parsecs of AGN. If this holds, nuclear 
star clusters grow only when SMBHs cannot. 

Finally, our model explicitly assumes that the nuclear 
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star clusters and bulges of dwarf ellipticals form in a quasi- 
spherical or at least a geometrically thick disc configuration 
of gas. If instead the gas is in a thin disc configuration before 
the onset of star formation, the feedback efficiency would be 
greatly reduced, and no significant bulge would be formed. 
Therefore our model does not apply to bulgeless spiral galax- 
ies. If central £ t ens of parsec of the se galaxies are fed via 
gaseous discs (e.g. jMilosavlievio!l2004l '). then the mass of the 
NCs need not saturate at the value given by equation [S] 

The assumption that star formation proceeds on a sin- 
gle dynamical timescale is a lower limit on the time actu- 
ally required. In fact it is more likely that star formation 
in the bulge takes several dynamical times to complete. In- 
deed, observationally we know that Giant Molecular Clouds 
in the Milky Way must be contracting much slower than 
dynamical collapse (|Zuckerman fc Palmed IT974T I to explain 
the low star formation efficiency in the Galaxy, pres umably 
due t o feedback by star formation inside the clouds (|McKeej 
1 19891 ). We would therefore expect a transition regime around 
a ~ 150 km s -1 where galaxies may contain either nuclear 
clusters or black holes. This boundary region extends over 
a factor of ~ 2 in a. In this region the competition between 
NCs and SMBH depends on the detail of gas deposition in 
the inner region of the galaxy and perhaps the merger his- 
tory of the galaxy. 

The picture we have presented is necessarily very sim- 
plified. One would like to include effects such as a realistic 
galaxy bulge potential, density inhomogeneities, and possi- 
ble cooling effects. Further, the changeover between NC and 
SMBH-dominated bulges depends on the merger history of 
the galaxy. For all these reasons a numerical treatment of 
this picture is desirable. The fact that our own Galaxy ap- 
pears to lie in the regime where the merger history may play 
an important role should make such studies rewarding. 
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